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In three-terminal NSN hybrid structures the influence of additional barriers on the nonlocal con- 
ductance and on current cross-correlations is studied within a scattering theory. In metallic systems 
with additional barriers and phase averaging, which simulate disordered regions, local processes can 
be enhanced by reflectionless tunneling but this mechanism has little influence on nonlocal processes 
and on current cross-correlations. Therefore Cooper pair splitting cannot be enhanced by reflec- 
tionless tunneling. On the contrary, in ballistic systems, additional barriers lead to Fabry-Perot 
resonances and allow to separate the different contributions to the conductance and to the current 
cross-correlations. In particular, crossed Andreev processes can be selectively enhanced by tuning 
the length or the chemical potential of the interbarrier region. 
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I. INTRODUCTION 



Transport in normal metal-superconductor-normal 
metal (N a SNh) hybrid three-terminal nanostructures of- 
fers a special interest, owing to the suggestion of pro- 
ducing split pairs of spin-entangled electrons from a 
superconductor^^. This is possible when the size of 
the region separating the N a S and the N^S interface 
becomes comparable to the superconducting coherence 
length. This circumstance indeed allows coherent pro- 
cesses involving two quasi-particles, each simultaneously 
crossing one of the two interfaces 3 - - — . Much effort has 
been devoted to the theoretical understanding and to 
the experimental observation of such a Cooper pair split- 
ting effect. In a transport experiment where electrons 
are difficult to measure one by one — contrarily to the 
similar production of entangled photons — , one should 
extract the relevant information from steady transport 
measurement, c. g. the current-voltage characteristics 
(conductance) and the cumulants of the current fluctu- 
ations (non-equilibrium current nois o 10 ' 11 and its count- 
ing statistics^). In practice, the conductance and the 
second-order cumulant (the shot noise and the current 
cross-correlations between the two current terminals N a , 
N(,) are the quantities to be extracted from experiments. 
Indeed, due to Fermi statistics, the "partition" noise cor- 
relations at a three-terminal crossing of normal metal 
contacts are negative^, manifesting the antibunching 
properties of individual electrons. If instead one contact 
is made superconducting, the cross-correlations may be- 
come positive, indicating the splitting of Cooper pairs^. 

There are two basic nonlocal processes occurring at the 
double NS interface: Crossed Andreev Reflection (CAR) 
which alone leads to Cooper pairs splitting into separated 
electrons bearing opposite spins (for a spin singlet super- 
conductor), and Elastic Cotunneling (EC) which alone 
leads to (spin-conserving) quasi-particle transmission be- 
tween the normal contacts, across the superconducting 
gap£. For tunnel contacts, at lowest order in the bar- 
rier transparencies, those two processes arc simply re- 
lated to the conductance, leading to positive (resp. nega- 



tive) conductance and current cross-correlations for CAR 
(resp. EC) processes. Indeed Bignon etalr^- showed 
that for tunnel contacts the linear dependence of the cur- 
rent cross-correlation on the voltages applied to the con- 
tacts N a , N;, allows to separately track the amplitudes 
of CAR and EC. Due to the expected compensation of 
the opposite CAR and EC conductance components at 
low transparencies, ferromagnetic contacts are required 
to detect CAR and EC from the conductance with tun- 
nel contact o 6 ' 7 ' 13 . Yet, such polarizations arc not easily 
achievable, moreover, if one is interested in producing 
spin-entangled electrons in a nonlocal singlet state, one 
should of course not spin polarize the contacts. 

Yet, as regards experiments, the situation for (ex- 
tended) tunnel barriers looks more complicated than 
given by a simple tunnel models. In addition, zero- 
frequency noise measurements involve weak signals and 
therefore they require more transparent interfaces, so 
that the arguments based on a tunnel model cannot be 
directly applied. For more transparent contacts, it has 
been found theoretically^^ - — that the nonlocal conduc- 
tance is dominated by EC thus is negative, which leaves 
the current cross-correlations as the only possible proof 
of Cooper pair splitting processes, provided one controls 
the voltages on both contacts. Contrarily to conductance 
measurements^, cross-correlations have led to few ex- 
perimental result o 20 ' 21 . At the theoretical level, the de- 
pendence of the cross-correlations on the contact trans- 
parency is not yet fully understood. Further analysis is 
required if one wishes to extract information from the 
sign and the amplitude of both the nonlocal conductance 
and the current cross-correlations in a system where the 
contacts are non-fcrromagnctic metals and the interface 
transparencies are rather large. 

In view of the current experiments on metallic struc- 
tures, the main question is therefore: Can the cross- 
correlations be positive, and if the answer is yes, is this 
a signature of Cooper pair splitting ? Previous work 
on a NSN structure^ showed that the cross-correlations 
can indeed be positive at large transparencies, although 
the nonlocal conductance is negative. The origin of this 
somewhat surprising result was not fully elucidated. Fur- 
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ther worls— showed that the sign of the cross-correlations 
indeed changes with the transparency of the interfaces, 
being positive at low transparency, negative at inter- 
mediate transparency and positive again at high trans- 
parency. While the positive sign at low transparency is 
clearly ascribed to Cooper pair splitting, it was shown 
that the positive sign at high transparency should not 
be interpreted in the same way. Indeed, at high trans- 
parency, CAR processes do not dominate either in the 
conductance or in the noise. Instead, the positive cross- 
correlations should be ascribed to higher order processes, 
in particular those involving local Andreev reflection 
(AR) on one side, and the opposite process on the other 
side, a process equivalent to exchanging a pair of elec- 
trons between the two normal contacts. 

These results show that transparent interfaces are 
detrimental to Cooper pair splitting. A possible alter- 
native strategy to boost Cooper pair splitting is to take 
advantage of multiple scattering to enhance the CAR 
process. Indeed, at a single NS interface, it was shown 
that disorder in the N region, or multiple scattering at 
a clean NN;S double interface, can strongly enhance An- 
dreev reflection, a mechanism nicknamed "reflectionless 
tunneling"— ~—. With a disordered N' region, it holds 
below the Thouless energy, where the mutual dephasing 
of electrons and Andreev-reflected holes is negligible. In 
the case of a clean double NN;S interface, maximum An- 
dreev transmission is obtained by balancing the trans- 
parencies T„„ and T ns of the NN; and N/S interfaces, 
such that the N; region acts as a resonant cavity. Melsen 
and Beenakker— performed an average on the modes in- 
side N; in order to mimick a disordered region. One can 
wonder whether a similar mechanism could enhance non- 
local Andreev reflection e.g. boost CAR compared to 
AR and EC. For this purpose we consider in this work 
a set-up N a N;SN r Nh with a quadruple interface. This 
system can be either considered as a modclization of a 
true one-dimensional hybrid structure with clean (bal- 
listic) normal and superconducting elements separated 
by tunnel barriers, for instance using carbon nanotubes 
or semiconducting nanowires. Or it can mimick a NSN 
system where the N metals are disordered close to the 
S interface. A scattering theory is performed in a one- 
dimensional geometry, varying the transparencies of the 
barriers and the width of the superconductor. 

One important result of this work is the following : 
If averaging independently the modes in the left and 
the right regions Nj, no boosting of the CAR process 
is obtained. Indeed, "reflectionless tunneling" is a quan- 
tum coherent process which demands that the Andrccv- 
reficcted hole retrace the path of the electron, by scat- 
tering on the same impurities. On the contrary, with 
nonlocal Andreev reflection, the electron and the trans- 
mitted hole sample different disorders and no coherence 
is obtained. Thus, the common reflectionless tunneling 
mechanism cannot enhance Cooper pair splitting in such 
a disordered one-dimensional geometry. Interestingly, we 
show that if on the contrary the N; regions can be tuned 
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Figure 1: Schematic of the model. 



at a Fabry-Perot resonance for the specific electron and 
Andreev-transmitted hole wavevectors, then CAR can be 
efficiently enhanced, moreover CAR and EC processes 
can be filtered at will by tuning the energy or the Fermi 
wavevector of both N/ region. The feasibility of such a 
proposal is discussed at the end. 

Section [TT| presents the model and section [TTTI the scat- 
tering theory of the N a N/SN r N{, system. Section |IV] pro- 
vides the results obtained by averaging over channels in 
N; and N r , in the spirit of Ref. Section [V] describes 
transport in the case where N; and N r are clean metallic 
"cavities". Section [VTl concludes by discussing a possible 
implementation of a controllable Cooper pair splitter. 



II. THE MODEL 

We study a one-dimensional model of a symmetri- 
cal three-terminal normal metal-superconductor-normal 
metal hybrid structure depicted in Fig. [T] The central 
superconducting electrode is grounded, the normal ter- 
minals can be biased with voltages V a and Vb- The length 
R of the superconducting electrode is comparable to the 
superconducting coherence length, which makes the ob- 
servation of nonlocal effects possible. The interfaces be- 
tween the normal metal and the superconducting elec- 
trodes are modeled by barriers with transparencies T/ ns 
and T rsn . In both normal metal electrodes there is an 
additional barrier at a distance Li (respectivly L r ) from 
the normal metal superconductor interface with trans- 
parency Ti nn (respectivly T rnn ). 

The system can be described by a 4 x 4 scattering 
matrix sff where Latin indices run over the normal 
electrodes a and b and Greek indices over electrons e and 
holes h. For energies lower than the superconducting 
gap, there are no single electron or hole states in the 
bulk of the superconducting electrode. Consequently, 
there is no scattering of electrons or holes in or out 
of the reservoir in the superconducting electrode and 
the latter does not appear in the scattering matrix. 
However, the scattering theory implicitly assumes that 
the superconductor is a reservoir of Cooper pairs. It 
is thus essential that the superconducting electrode is 
grounded. The transformation of quasi-particles into 
Cooper pairs is taken into account by the correlation 
length £ which sets the scale of the damping of the 
electron and hole wavefunctions in the superconductor. 
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The elements of the scattering matrix are evaluated 
from the BTK approach 2 ^ (see Appendix IA"]) . As the en- 
ergy dependence of the calculated quantities is relevant, 
we do not use the Andreev approximation valid in the 
limit of zero energy where the electron and hole wavevec- 
tors are set to the Fermi wavelength, but instead keep the 
full expressions for the wave vectors. The calculation of 
the average current and current cross-correlations rely on 
the formulas obtained by Anantram and Datta in Ref. H3: 



fc£{a,b} a,P<E{e,h} 



x J dE 

2e 2 

Sij = — H sgn(a)sgn(/3) 

k,l£{a,b} a,P,-f,5£{e,h} 
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x dE A kl; is(ia,E)Ai SM (ia,E)f k ^(E)[l- fi S (E)} 



fects are studied in section [V] In a higher-dimensional 
system with more than one transmission mode, the os- 
cillations in the different modes arc independent and are 
thus averaged out. Multidimensional behavior can be 
simulated qualitatively with a one-dimensional system by 
averaging all quantities over one oscillation period: 



X(Li,R, L r ) 



(3) 
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drX(k,r, l r ). 



This procedure is appropriate to describe metallic sys- 
tems. These averaged quantities are studied in section 



III. COMPONENTS OF THE DIFFERENTIAL 
CONDUCTANCE AND THE DIFFERENTIAL 
CURRENT CROSS-CORRELATIONS 



with A kl js{ia,E) = 8 ik 8 a 5 ai 8 a s - s" fe s-f, 
sgn(a = e) = 1 sgn(a = h) = —1, fj a the oc- 
cupancy factors for the electron and hole states 
in electrode i, given by the Fermi function where 
the chemical potential are the applied voltages 
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In this one-dimensional model, both current and noise 
are highly sensitive to the distances Li, R, L r between 
the barriers: They oscillate function of these dis- 
tances with a period equal to the Fermi wavelength Xp <C 
Li,R,L r . In a ballistic system, the multiple barriers act 
like a Fabry-Perot interferometer and the interference ef- 



An electron, arriving from one of the normal metal 
reservoirs at the interface to the superconductor, can 
have four different destinies: It can be reflected as 
an electron (normal reflection (NR)), or reflected as a 
hole (Andreev reflection (AR)), or transmitted as an 
electron (elastic cotunneling (EC)) or transmitted as a 
hole (crossed Andreev reflection (CAR)), and similarly 
for holes. The corresponding elements of the scattering 
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The current in electrode N a given by Eq. (TTJ) can natu- 
rally be divided into AR, CAR and EC contributions (to 
calculate this expression the unitarity of the scattering 
matrix has been used): 



I a = ^ J dE [(\sf a (E)\ 2 + \s h a Z(E)f) (f ae (E) - f ah (E» 

local Andreev reflection 

+ \.s:t(E)\ 2 (f ae (E) - f he {E)) + \s h a £(E)\ 2 (f bh (E) - f ah (E)) 

clastic cotunneling 

+ \sf b {E)\\f ae {E) f bh {E)) + \s%(E)\ 2 (f be (E) f ah {E))} . (4) 

crossed Andreev reflection 

In the following, we focus on i) the differential conductance in the symmetrical case where V a = V& = V and the 
current I a is differentiated with respect to V, and ii) the differential nonlocal conductance in the asymmetrical case 
where V a = and the current I a is differentiated with respect to V b - In the zero temperature limit, only the nonlocal 
processes, CAR and EC, contribute to the nonlocal conductance: 

= ~ [Kt(\e\V b )\ 2 + I^He|Vi)| 2 ] + j [| s ^(-|e|V b )| 2 + |^(|e|F 6 )| 2 ] , (5) 

clastic cotunneling crossed Andreev reflection 
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while the symmetric case contains local Andrccv reflection and crossed Andreev reflection: 

= T [(Ka(NV0| 2 + l^(M^)| 2 ) + (l^(-|e|^)| 2 + |*£(-|e|V)| 2 )] 
v a =v b =v J± v , 

local Andreev reflection 

+ j [H1(\e\V)\ 2 + \s h a t(\e\V)\ 2 ) + (K h b (~\e\V)\ 2 + I^H^)| 2 )] (6) 

V v ' 

crossed Andrccv reflection 

Let us now perform a similar analysis for the current cross-correlations. We study only the zero temperature limit, 
where fk^(E)[l — fis{E)] is zero if k = I and 7 = 8 and the current cross-correlations are: 

S ab (T = 0) = ^ £ sgn(a)sgn(/3) f dEs^ s<*J sft f kl (E)[l f ls (E)} (7) 

k,le{a,b} a,P,-Y,8G{e,h} 



dla 

dV 



Every summand in S a b contains the product of four ele- 
ments of the scattering matrix. As pointed out in Refs. 
[28ll29l . in difference to the situation for the current, it is 
impossible to combine those matrix elements to absolute 
squares. Let us now try to classify the contributions of 
the noise as we did above for the current. We find that 
no summand consists of only one kind of elements of the 
scattering matrix. Every element consists of two local 
elements (NR or AR ) and two nonlocal elements (CAR 
or EC). Either the two local elements and the two non- 
local elements are identical, that gives the components 
EC-NR, CAR-NR, EC-AR, CAR-AR, or all four matrix 
elements belong to different categories and we will call 
these summands MIXED. Sometimes, it is useful to di- 
vide MIXED further as a function of its voltage depen- 
dence (see Appendix [Bj . As the formulas for the current 
cross correlations are lengthy, they are relegated into ap- 
pendix [B] 

For the interpretation of current cross-correlations, the 
global sign plays an important role. For positive applied 
bias voltages, the differential current cross-correlations 
carry the same sign as the current cross-correlations. For 
negative applied voltages current, cross-correlations and 
differential current cross-correlations have opposite signs. 
To avoid confusion, we only show pictures of the dif- 
ferential current cross-correlations calculated for positive 
bias voltages (and thus negative energies E = —\e\V). 
Due to the electron-hole symmetry of the model, differ- 
ential current cross-correlations calculated for negative 
bias voltages are up to a global sign identical to the ones 
calculated at positive bias voltage. For small bias volt- 
ages, current cross-correlations depend linearly on volt- 
age. Thus, current cross-correlations and differential cur- 
rent cross-correlations show the same qualitative behav- 
ior if studied as a function of the interface transparency 
or the distance between barriers. 

In Ref. Il8l . we performed a similar analysis of cur- 
rent cross-correlations in terms of Green's functions. For 
the relations between these two classifications see Ap- 
pendix [Cj 



Bignon et al^ have studied current cross-correlations in 
the tunneling limit. They find that noise measurements 
in the tunneling limit can give access to the CAR and 
EC contribution of the current. We have just seen, that 
at least two processes are involved in every component 
of noise, but the contributions of noise they calculate fall 
into the categories EC-NR and CAR-NR. In the tunnel- 
ing limit the NR contribution is so close to one, that it 
can be neglected and what remains is very similar to the 
current contributions. 



IV. THE AVERAGED SYSTEM 
A. Positive Cross-Correlations without CAR 

In fermionic systems without superconducting elec- 
trodes, current cross-correlations arc, as a consequence 
of the Fermi-Dirac statistics, always negativ o 10 i 27 . In 
superconductors, on the contrary, the coupling of elec- 
trons into Cooper pairs makes positive cross-correlations 
possible^. If a CAR process is interpreted as the splitting 
of a Cooper pair into two electrons leaving the supercon- 
ductor in different electrodes, positive cross-correlations 
are its logical consequence. However, the CAR process 
is not the only one which can lead to positive cross- 
correlations. Let us investigate in more detail the in- 
fluence of the different processes on the current cross- 
correlations in a NSN-system, i. e. in a system without 
the additional barriers in the normal conducting elec- 
trodes. 

The black line in Figure [D shows the averaged dif- 
ferential current cross-correlations for symmetric bias 
(V = V a — Vb). The total current cross-correlations have 
already been published in Ref. [H, but here, Figure [5] 
shows in addition the different parts which contribute 
to the total current cross-correlations. The total cross- 
correlations are positive for high interface transparencies 
and for low interface transparencies. As we have already 
argued in Ref. [IH, the positive cross-correlations at high 
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Figure 2: Averaged differential current cross-correlations for 
a symmetrical biased (V = V a = Vt -C A/|e|) NSN-system as 
a function of the transparency of the interfaces T; ns = T rm . 
The positive cross-correlations at high interface transparency 
are due to the EC-AR process, represented by a dotted line. 



interface transparencies are not due to CAR. Indeed, at 
very high transparencies, only processes which conserve 
momentum can occur, since there are no barriers which 
can absorb momentum. CAR processes do not conserve 
momentum: if e.g. an electron arrives from the left hand 
side carrying momentum the hole that leaves at the 
right hand side carries momentum — kp. As the cross- 
correlations do not tend to zero even for very high trans- 
parencies, they cannot be due do CAR. Indeed, if we plot 
the different components of the noise introduced in the 
last section separately, we see that positive current cross- 
correlations at low interface transparencies are a conse- 
quence of a large CAR-NR component and therefore a 
consequence of CAR processes. But the positive current 
cross-correlations at high interface transparencies have a 
different origin: a large positive EC-AR contribution. 

We can put the contributions to the current cross- 
correlations into two categories with respect to their sign, 
which is independent of the interface transparency. EC- 
NR, CAR-AR, MIXED2 and MIXED4 carry a negative 
sign, CAR-NR, EC-AR, MIXED 1 and MIXED3 carry a 
positive sign. The current can either be carried by elec- 
trons I e or by holes J . The sign of the different con- 
tributions to the current cross-correlations depends on 
whether only currents of the same carrier type are corre- 
lated^ ((A^A4 e ) + (Ai£A^ l )+a <-> b), which is the case 
for EC-NR, CAR-AR, MIXED2 and MIXED4 and leads 
to a negative sign; or whether electron currents are corre- 
lated with hole currents ((Ai^AI^) + (AI^Ai^) +a^b), 
which is the case for CAR-NR, EC-AR, MIXED 1 and 
MIXED3 and leads to a positive sign. In purely nor- 
mal conducting systems, the electron and hole currents 
are uncorrelated, only correlations of the same carrier 
type contribute to the current cross-correlation and lead 
to a negative sign. The sign of the total current cross- 
correlations is a consequence of the relative strength of 
the different parts of the current cross-correlations, which 
depends on the interface transparency. 



B. Multiple Barriers 

In the last paragraph, we learnt that positive cross- 
correlations due to CAR can only be found in the tun- 
neling regime. But in the tunnel regime the signals are 
quite weak. The conductance over an NS-tunnel junction 
can be larger for a "dirty" normal conductor contain- 
ing a large number of non-magnetic impurities, where 
transport is diffusive, than for a clean normal conduc- 
tor, where transport is more ballistic, by an effect called 
rcflcctionlcss tunnelin g 23 i 24 . It is thus natural to ask if 
a similar effect could also enhance conductance and cur- 
rent cross-correlations in a three terminal NSN-structure. 
To answer this question, we use the model of Melsen 
and Beenakker— where the disordered normal conduc- 
tor is replaced by a normal conductor with an additional 
tunnel barrier leading to an NNS structure. Duhot and 
Melir>i£ have studied the influence of additional barri- 
ers on the nonlocal conductance in three terminal NSN- 
structures. They find that two symmetric additional bar- 
riers enhance the nonlocal conductance. 

First, let us get a deeper understanding of their re- 
sult by calculating the AR, CAR and EC components 
of the current separately. Afterwards, we will study 
the influence of additional barriers on the current cross- 
correlations. 

Figure [3] shows the averaged conductivity in the sym- 
metrical bias situation V a = Vb <C A/|e| and in the 
asymmetrical voltage case V a — 0, Vb <C A/\e\ for a 
superconducting electrode much shorter than the coher- 
ence length {R = 0.250- The sum of the AR, CAR 
and EC components, traced in black, features in both 
cases an extremum. Yet, looking at the behavior of there 
components, we see that they arise from different mech- 
anisms. Let us first have a look at the symmetrically 
biased case. Without the additional barriers, i.e. in the 



limit Ti 



T 



1, the contributions of AR and 



CAR are similar in magnitude. The EC component is 
completely suppressed, since it is proportional to the dif- 
ference of the applied voltages. The introduction of two 
additional barriers increases the AR component by about 
a factor 30. The shape of the curves is similar to the 
one of the NNS structure derived analytically by Melsen 
and Beenakker—, which can be recovered exactly by in- 
creasing the length of the superconducting electrode far 
beyond the coherence length. However, the CAR curve 
stays almost constant over a wide range of values of bar- 
rier strength of the additional barriers, and it eventually 
vanishes when the transparencies go to zero. 

In the asymmetrical voltage case V a = 0, the AR com- 
ponent is zero as it is proportional to the local voltage 
V a - Like in the first case, the additional barriers have lit- 
tle influence on the CAR component, except for the fact 
that it tends to zero for vanishing transparency. Over 
a wide range of barrier strength values, the EC compo- 
nent is identical in amplitude, but opposite in sign to the 



CAR component. For small T;, 



T r nn values the EC 



component displays a small extremum, but it is much 
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Figure 3: Averaged differential conductance in the limit of 
zero energy in a) the symmetrical bias situation V a = Vb -C 
A/|e| and b) in the asymmetrical case V a = 0, Vb <C A/|e[ for 
a superconducting electrode much shorter than the coherence 
length (R = 0.25£) as a function of the transparencies of 
the additional barriers Ti n „ = T rnn . The barriers next to the 
superconductor are in the tunnel regime (Ti ns = T rS n = 0.01). 



less pronounced than the maximum of the AR compo- 
nent of the first case. In the limit Ti nn = T rnn — > 1, 
the EC component tends more slowly to zero than the 
CAR component which as a consequence has a minimum 
in the total conductance. The fact that the conductance 
maxima in the symmetrical bias case and in the asym- 
metrical bias case have different origins can also be illus- 
trated by studying their energy dependence, depicted in 
Figure 2J The enhancement of the AR component of the 
conductance in the symmetrical biased case disappears 
completely with increasing bias voltage. The extremum 
of the conductance in the asymmetrical biased case de- 
creases slightly with increasing bias voltage, but only up 
to a certain voltage value, then it saturates. 

Why is the AR component enhanced by the additional 
barriers, but not the EC or CAR components? Refiec- 
tionless tunneling is believed to occur because the elec- 
trons and holes get localized between the double barrier 
and have therefore a higher probability to enter the su- 
perconductor. There is a double barrier on the right and 
on the left hand side, so this localization should also hap- 
pen to the particles involved in EC or CAR processes. 
The answer is, that in the AR cases the incoming elec- 
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Figure 4: Total averaged differential conductance, as in Fig- 
ure [3] but at different energies. 



tron and the leaving hole see the same environment. In 
the EC and CAR process, the incoming particle sees the 
environment on one side of the superconductor and the 
leaving particle the environment on the other side. The 
energy dependence of the conductance enhancement of 
the AR component demonstrates that reflectionless tun- 
neling can only occur if the outgoing hole can trace back 
the way of the incoming electron. At low bias voltages 
where the electron and the hole have nearly the same 
energy and wavevector, reflectionless tunneling occurs. 
At higher bias voltage, electrons and holes have different 
wavevectors and the reflectionless tunneling peak disap- 
pears. The integrals over the phases between the addi- 
tional barriers on the left- and on the right-hand side 
have, of course, been taken independently. There is no 
reason to think that the channel mixing, which is emu- 
lated by the integrals, on the left- and on the right-hand 
side are coupled. To verify this scenario, let us couple 
the two integrals in an gcdankcncxpcriment. We set the 
distance Li between the two left-hand side barriers to be 
equal to the distance L r between the two right-hand side 
barriers and perform only one integral over L = L\ = L r . 
The result is shown in Figure Now, the CAR and the 
EC component are also enhanced. Still, the effect on the 
EC component is larger and EC dominates the nonlocal 
conductance. 
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Figure 5: Gedankenexperiment with coupled integrals Ti ns — 
Trsn = 0.01, R = 0.25£: Now, also the EC and the CAR 
component are enhanced by reflectionless tunneling. 



Let us turn back to independent averaging and have a 
look at the current cross-correlations shown in figure [6] 
In the symmetrical bias case, the additional barriers do 
not lead to an enhancement of the signal. The noise 
is dominated by the CAR-NR component, and, as we 
have seen above, CAR is not influenced by reflectionless 
tunneling. The EC-AR component is amplified by the 
additional barriers, because the AR amplitude describing 
a local process is amplified. This leads to a small shoulder 
in the total cross-correlations. But since we are in the 
tunnel regime and the leading order of CAR-NR is T 2 
while the leading order of EC-AR is T 4 , the influence 
of the EC-AR-component is to small to lead to a global 
maximum. 

In the asymmetrical bias case V a = 0,Vj, <C A/|e|, on 
the other hand, the additional barriers enhance the sig- 
nal. But the cross-correlations are dominated by EC-NR 
and are therefore negative. We conclude that in a phase- 
averaged system, additional barriers only enhance the 
AR-componcnt, a local processes. It cannot help to am- 
plify nonlocal signals. Thus, positive cross-correlations, 
when due to nonlocal processes cannot be enhanced with 
additional barriers if an average over the length has to 
be done. 



V. THE BALLISTIC SYSTEM 

In a last step, let us examine if the additional barriers 
can enhance the current cross-correlations in a ballistic 
one-channel device, where no integrals over the length of 
the normal conducting channels need to be done. We 
nevertheless maintain he average over the central su- 
perconducting electrode. One can indeed imagine that 
the normal conducing channels are formed by single- wall 
nano-tubes, where transport is ballistic, while the super- 
conducting central electrode is made from more disor- 
dered aluminum. 
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Figure 6: Averaged differential current cross-correlations in 
a) the symmetrical bias situation V a = H <C A/|e| and b) 
the asymmetrical bias case V a = 0, Vb <S A./|e| for a su- 
perconducting electrode shorter than the coherence length 
(R — 0.25£) as a function of the transparencies of the ad- 
ditional barriers T; nn — T rnn . The barriers next to the super- 
conductor are in the tunnel regime (Ti ns = T rsn = 0.01). 



A. Resonances in a Fabry-Perot Interferometer 

Assuming ballistic transport in the normal conducting 
parts, the NNSNN-system consists in principle of two 
Fabry-Perot interferometers, one at the left- and one at 
the right-hand side of the superconductor. For the final 
result the entire system has to be taken into account. 
But thinking of the two Fabry-Perot interferometers in- 
dependently already gives a fair idea of where to expect 
resonance peaks, and what influence they have on the dif- 
ferent components of the conductance and of the current 
cross-correlations. 

In an optical Fabry-Perot interferometer consisting of 
two parallel mirrors, the transmitted intensity results 
from the interference of light which has undergone mul- 
tiple reflections between the two mirrors. The transmis- 
sion of a Fabry-Perot interferometer is a function of the 
phase difference between a beam which has been n-times 
and a beam which has been (n + l)-times reflected at 
both interfaces. The same is true for a Fabry-Perot in- 
terferometer for electrons. An electron going once back 
and forth between the two barriers at the left-hand side 
of the NNSNN-system acquires a phase </> e = 2q + Li, 
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Figure 7: Schematics of the contributions to the conductance 
in the symmetrically biased case. The filled states are de- 
picted in the three reservoirs, together with the electron and 
hole currents between the barriers, a) AR contribution for 
positive bias, b) and c) The two possible CAR contributions 
for positive bias, d) Summary of the processes for negative 
bias. If negative voltages are applied, the roles of electrons 
and holes are interchanged, but the resulting conductance 
stays unchanged. 



a hole acquires a phase <f> h = 2q~Li, where q + is the 
wavevector of the electrons and q~ the wavevector of the 
holes. For low voltages, where q + ~ q~ , we expect elec- 
tron and hole resonances to happen for the same res- 
onator len gths. I f the voltage is increased, the wavevec- 
tor q + = \ / 2mh 2 y/Jj, — \e\V decreases and the wave vec- 
tor q~ = ^2mh 2 ^/ /it + \e\V increases, shifting the reso- 
nance for electrons to higher and the resonance for holes 
to lower values of L;. The width of the resonances of a 
Fabry-Perot interferometer decreases as the reflectivity 
of the interfaces increases. 



1. Symmetrically applied bias 

For a symmetrically applied bias, the differential con- 
ductance consists of an AR and a CAR component. If 
a positive bias is applied to the normal electrode N a , 
the current I a contains an AR contribution, where a hole 
enters and an electron leaves the superconductor (see fig- 
ure [7^). At low energies, the AR contribution will have a 
maximum for values of L\ where electrons and holes are 
both in resonance. At higher energies, the resonances 
for electrons and for holes occur at different values of Li . 
If the barriers are in the tunnel regime, the resonance 
peaks are sharp. The resonance peaks for electrons and 
holes stop overlapping and the AR contribution disap- 
pears. If the interface transparency is increased, the res- 
onance peaks broaden. Even if the electron and the hole 
wavevectors are significantly different, there is an over- 
lap between corresponding resonances leading to a double 
peak in the AR contribution. AR is a local process, and 



a) 
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Figure 8: Schematics of the contributions to the conductance 
in the asymmetrically biased configuration. The filled states 
are depicted in the three reservoirs, together with the current 
of electrons and holes between the barriers, a) EC contribu- 
tion for positive bias, b) CAR contribution for positive bias. 



we expect the AR contribution to the current I a to be 
mainly independent of the length L r . 

There are two possibilities for CAR process if the elec- 
trodes N a and Nb are symmetrically and positively bi- 
ased: Either a hole enters the superconductor at the left- 
hand side and an electron leaves the superconductor at 
the right-hand side (sec figure EJd); or the process occurs 
the other way around (see figure [Tj:), the hole enters at 
the right-hand side and the electron leaves at the left- 
hand side. In a CAR process, both sides are involved 
and the position of the resonances will be a function of 
both Li and L r . At low energies, the resonances for the 
two CAR processes will be superimposed and occur at 
length values Li = L r mod n/kp- In a CAR process, 
the incoming hole and the outgoing electron are spatially 
separated. It is therefore possible to reach the optimal 
resonance condition for electrons and for holes even if 
the two kinds of carriers have different wavevectors. For 
higher bias values, the conductance maximum due to the 
CAR process sketched in figure[7}3 will be shifted to lower 
values of Li and higher values of L r . It will be the other 
way around for the conductance maximum due to the 
second CAR process. Therefore, we expect the resonance 
peak to split up into two peaks. 

CAR and AR both have a positive contribution to the 
differential conductance. This can also be seen in figure[7] 
The current I a is defined as positive if positive charges 
enter the superconductor. On the left-hand side, all elec- 
trons leave the superconductor and all holes enter it. If 
a negative bias is applied instead of a positive one, elec- 
trons and holes will interchange their roles (sec figurc[7Jl). 
The resonance positions and the sign of the conductance, 
however, will not change because q + (V) = q~(~V). 



2. Asymmetrically applied bias V a — 0, Vb > 

Let us now focus on the bias configuration V a = 0, Vt > 
used to measure a nonlocal conductance. The latter 
consists of an EC and a CAR contribution. EC processes 
involve the same type of charge carriers on both sides of 
the superconductor. Therefore, the resonances for EC 
will always occur for diagonal values L/ = L r mod w/kF- 
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At positive bias, EC is carried by holes (see figure^) and 
if the voltage is increased, the resonances in Li and L r 
are simultaneously shifted to lower values. EC leads to a 
negative differential conductance as the holes move out 
of the superconductor and the current is defined positive 
when entering the superconductor. 

From the two kinds of possible CAR processes in the 
symmetrically biased case, only the second one con- 
tributes to the nonlocal conductance. At low energies the 
EC and the CAR peak will be superimposed. At higher 
energies, the CAR peak moves away from the diagonal. 
The application of a negative bias again interchanges the 
roles of electrons and holes, but it does not influence the 
resonance positions. The only difference between a posi- 
tive and a negative applied bias is the global sign of the 
current cross-correlations. 

If one wishes to use a Fabry-Perot interferometer as a 
filter for the different processes contributing to the con- 
ductance and to the current cross-correlations, we expect 
that it will be essential to operate the system at energies 
where the wavevectors of electrons and holes are signif- 
icantly different. In the next section, numerical results 
will be analyzed, which will turn out to be in qualitative 
agreement with the scenario discussed in this section. 
We start with the tunnel regime, then investigate how 
the contrast diminishes with increasing interface trans- 
parency and eventually have a closer look at a system 
with intermediate interface transparencies. 



B. The Tunnel Regime 

Let us now discuss the calculated conductance and the 
current cross-correlations. In a first step, the barriers 
between the superconducting and the normal regions are 
chosen to be in the tunnel regime (T/ ns = T rsn m 0.01) 
in order to achieve sharp resonances. The distances 
Li and L r are of the order of ten times the coherence 
length £. This distance allows to observe phase differ- 
ences Acj> = \4> e — (j) h \ of the order of 2tt between elec- 
trons and holes already at energies much smaller than 
the superconducting gap. 

As for the averaged system, we will study four dif- 
ferent quantities: the conductance in the symmetrically 
biased case, the nonlocal conductance, where the nor- 
mal electrode in which the current is calculated is at the 
same potential as the superconductor, and the current 
cross-correlations in the same two bias configurations. 
The results are presented in the figures |H [TU1 [TT] and [T2] 
The panels a) and c) always show a three-dimensional 
plot of the conductance (respectively the current cross- 
correlations) as a function of the distances between the 
barriers Li and L r . The three-dimensional plot in panel 
a) always depicts the situation in the limit of zero en- 
ergy. The one in panel c) refers to higher energies. The 
three-dimensional plots are completed by cuts through 
the resonances. These cuts show the different contribu- 
tions in different colors and have a higher resolution. 




Figure 9: Differential conductance for symmetrically applied 
bias V = V a = Vi in the tunneling limit Ti ns = T rsn = 
0.01, Ti nn — T rnn = 0.26, a) as a function of L r and Li 
at low voltage \e\V = 1 ■ 10~ 4 A, b) cut through a) along 
one of the "mountain ranges" (at Li = 3 183.97r/fcF) at higher 
resolution, c) as a function of L r and Li at higher voltage 
\e\V = 6.25 ■ 10" 2 A, d) cut through c) at U = 3184.1tt/A:f at 
higher resolution. 




b) 



Figure 10: Differential nonlocal conductance (V a = 0) in the 
tunneling limit Ti ns = T rs „ = 0.01, Ti nn = T r „„ = 0.26, a) as 
a function of L r and Li at low voltage \e\Vb = 1 ■ 10 _4 A, b) 
cut through a) at L r = 3 183.97r/fcF at higher resolution, c) as 
a function of L r and L; at higher voltage \e\Vb = 6.25- 10 -2 A, 
d) cut through c) at L r = 3 183.67r/fcF at higher resolution. 



1. Conductance for a symmetrically biased system 

Figure [9^ shows the differential conductance 
dI a /dV\ v _y = y b>0 in the symmetrically biased case 
V a — Vb as a function of the distances L/ and L r in 
the limit V — > 0. There arc resonance "mountain 
ranges" for two values of Li and there are four peaks 
which stand out of the mountain ranges. The mountain 
ranges correspond to a L r -independent process and will 
therefore be due to AR. This is confirmed by figure [SJa 
showing in a cut along one of the mountain ranges the 
AR and the CAR component separately. Figure [SJd also 
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Figure 11: Differential current cross-correlations for V a = Vb 
in the tunneling limit T; ns = T rm = 0.01, T; nn = T„ m = 
0.26, a) as a function of L r and Li at low voltage \e\V = 
1 ■ 10~ 4 A, b) cut through a) at Li — 3183.97r/fcF at higher 
resolution, c) as a function of L r and Li at higher voltage 
\e\V = 6.25 ■ 10" 2 A, d) cut through c) at U = 3184.l7r/fc F at 
higher resolution. 
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Figure 12: Differential current cross-correlations for V a = 
in the tunneling limit Ti ns = T rsn = 0.01, T in „ = T r „„ = 
0.26. The three-dimensional plots for the current cross- 
correlations are quasi-identical to the ones for the nonlocal 
conductance in figure \Wl therefore only the cuts are shown, 
a) L r = 3183.97r/fc F , \e\V = 1 ■ 10~ 4 A, b) L r = 3183.67r/fc F , 
\e\Vb = 6.25 • HT 2 A 



shows that the mountain peaks are due to CAR. The 
relative distance of two subsequent resonance peaks in 
units of kp is tt in L r and in L\ direction. But the 
absolute resonance positions are not given by integer 
multiples of 7r as it would be the case in a simple 
optical Fabry-Perot interferometer with two identical 
mirrors. Figure [St shows the differential conductance 
at higher voltage (but still much smaller than the gap). 
As expected from section IV A[ the AR mountain ranges 
disappear because the electron and hole resonances do 
not overlap anymore. Every CAR peak splits up into 
two peaks corresponding to the two CAR contributions. 

The observation of the disappearance of AR and the 
splitting of the CAR peaks is already interesting in itself. 
As the peaks are uniquely due to CAR processes, the 
desired filtering of the different processes contributing to 
the conductance is achieved. A NNSNN-system operated 
with symmetrical bias at energies where AR is suppressed 
and tuned to a CAR resonance peak could be a good 
source of entangled electrons. 



2. Non-local conductance for V a = 0, Vb > 

Figure [TU] depicts the calculated differential nonlocal 
conductance for the biases V a = 0, V& > 0. Panel a) 
shows the nonlocal differential conductance dI a /dVh\ v _ 
as a function of the distances Li and L r at very low volt- 
age Vb- There are two resonance peaks on the diago- 
nal, where Li = L r , and two additional resonance peaks 
where L; = L r ± ir/kp. From the considerations in the 
last section, we expect the EC and the CAR resonances 
to appear at the same positions. As the peaks are nega- 
tive, the EC contribution is apparently stronger than the 
CAR contribution. Figure [TUb showing a cut through 
figure [TTJb at constant L;, where the total current is de- 
composed into its EC and CAR component, confirms that 
the peaks in the three-dimensional plot are indeed a su- 
perposition of EC and CAR. 

When Vb is increased, the EC and the CAR resonance 
move apart (see figure ITUb) . As expected, the EC peaks 
stay on the diagonals and at Li = L r ± ir/kF, but their 
positions are shifted to lower length values. The positive 
conductance peaks due to CAR appear at the same L r 
value as the EC resonance peaks, as both processes are 
carried by holes on the right hand side. On the left hand 
side CAR is carried by electrons whose wave vector q + 
decreases with increasing voltage V>. Consequently, the 
Li values of the peak position increase with increasing 
voltage. Again, the desired filtering of different processes 
is obtained. This voltage configuration is interesting as 
the conductance switches between positive and negative 
values if L r is changed. 



3. Differential current cross-correlations 

The differential current cross-correlations 
dS a b/dV\ v _y _ v depicted in figure QT] and 
dS a b/dVb\ v =Q depicted in figure [T^] behave very 
similarly to the corresponding differential conductance 
if one compares the CAR-NR component of the current 
cross-correlations with the CAR component of the 
conductance and the EC-NR component with the EC 
component. The only striking difference is the absence 
of the "mountain range" in the noise at low symmetrical 
bias. The similarity between current cross-correlations 
and conductance has already been discussed at the 
end of section IIIII in connection to the work of Bignon 
etalJ^: As the barriers are opaque, the most likely 
process to occur is normal reflection. Therefore, the 
elements of the scattering matrix corresponding to 
normal reflection are close to one and EC-NR is similar 
to EC and CAR-NR is similar to CAR. Only the AR 
component of the conductance has no counterpart in the 
current cross-correlations, because all elements of the 
current cross-correlations are composed of a local and a 
nonlocal process. An AR-NR component does not exist 
and the AR mountain range is absent in the noise. 
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Figure 13: Non-local conductance at \e\Vb = 6.25- 10 2 A and 
V a = for T in „ = Ti n3 = T r3n = T rnn = T as indicated in the 
figure, for a very short superconducting electrode R — 0.25£, 
at constant L r value, chosen to be in resonance so that the 
conductance in function of Li features CAR and EC peaks. 
The CAR peak splits up in two peaks. For T = 0.5 and larger, 
the CAR peak has completely disappeared. 



C. Increasing the Interface Transparency 



Figure 14: Non-local conductance at eVb = 6.25 ■ 10 2 A and 
V a — for T lnn = Ti ns = Trsn = T rnn — T as indicated 
in the figure, for a slightly longer superconducting electrode 
R = £, at constant L r value, chosen to be in resonance so 
that the conductance in function of Li features CAR and EC 
peaks. The EC and CAR peaks are sharper than in the case 
of the extremely short superconducting electrode depicted in 
figure [13] and there is no splitting. Even for barrier trans- 
parencies as high as T — 0.7 the signature of the EC and the 
CAR contribution is clearly visible. 



In the last section, we have studied the ballistic sys- 
tem in the tunnel regime where all resonance peaks are 
sharp. However, the systems in which electronic Fabry- 
Perot interferences have already been observed, for ex- 
ampl o 30 i 31 , were in the opposite regime of high inter- 
face transparency. For increased barrier transparencies, 
one expects the resonance peaks to become broader and 
the separation of the different peaks to become less pro- 
nounced. This section discusses whether it is necessary to 
stay in the tunnel regime or if the filtering effect survives 
to higher interface transparencies. 

To investigate systematically what happens when the 
interface transparency is increased, the nonlocal conduc- 
tance at eVh = 6.25 • 10 _2 A and V a = was calculated 
for different barrier transparencies, at constant L r value, 
which was chosen to be in resonance so that the con- 
ductance in function of Li would feature CAR and EC 
peaks. Figure [T31 shows this conductance for a short su- 
perconducting electrode R = 0.25£ and for different val- 
ues of Tl„n = Ti„s = T rs „ = T rnn = T. The CAR 
peak shows an interesting feature: It splits up into two 
peaks, while the EC peak does not. One possible ex- 
planation is the interplay between the resonance in the 
NNS double barrier and the resonance in the central su- 
perconductor, both involving Andreev reflections. The 
first one implies De Gennes-Saint- James bound states 3 -^, 



broadened by the presence of infinite reservoirs. The sec- 
ond are the Tomasch-Rowcll-McMillan^ 3 - resonances, due 
to a scattering between quasiparticle states with differ- 
ent wavevectors k + , k~ in S. In a NS bilayer, these two 
processes can altogether lead to a doubling of the reso- 
nance peaks^l. In the present case, such a doubling is 
observed in the CAR channel, but is absent in the EC 
channel where all wavevectors are approximately equal. 
Moreover, higher transmissions favour the EC channel, 
and the CAR peak gets weakened altogether. 

The amplitude of the quasiparticle wavefunction de- 
creases exponentially with the superconducting length. 
Therefore, one can hope to decrease the coupling between 
the two resonances by elongating the superconducting 
electrode. This actually works, as the plots in figure Q3] 
where R — £ show: The EC and CAR peaks are sharper 
than in the case of the extremely short superconduct- 
ing electrode and there is no splitting. Even for barrier 
transparencies as high as T = 0.7, the signature of the 
EC and the CAR contribution is clearly visible. Against 
the general trend, we have here a situation where a longer 
distance in the superconductor leads to a larger signal. 

In conclusion, it is not necessary to operate the 
NNSNN device in the tunnel regime in order to filter 
the different processes. Surprisingly, CAR resonances 
are better resolved if the distance in the superconducting 
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electrode is not too short. 

Let us now study a system for intermediate interface 
transparencies T = T/„„ = Ti„ s = T rsn = T rnn = 0.5 and 
R = £, which is clearly outside the tunnel regime but has 
still a good resolution between the different processes, in 
more detail. We proceed in the same way as in the tun- 
nel regime and analyze the conductance and the current 
cross-correlations in two different voltage configurations. 



1. Conductance for a symmetrically biased system 

Figure [TBI shows the differential conductance for a sym- 
metrically biased system. At low energies (figure ITSa ) 
the predominant structures are again the AR "mountain 
ranges" . But in contrast to the situation in the tunnel 
regime, there are no CAR peaks on top of the mountain 
ranges. Instead, at the positions we would have expected 
them, there are depressions. The cut along the moun- 
tain range (figure 115b ) shows that the AR component 
features a deep depression where L r = Lj. This depres- 
sion already exists in the tunnel regime (figure |9Jd) , but 
is much less pronounced there. The CAR component in- 
creases when L r approaches L; as if a peak was to form, 
but falls back to zero at L r = L\. We can understand 
these dips in the conductance from the competition of the 
different processes at the level of the scattering matrix. 
At L r = Li = 3 1 83.77t//cf, EC processes are in reso- 
nance (see section |VC2[) . As the probabilities for the 
four possible processes have to sum up to one, a large 
EC scattering matrix element is only possible at the cost 
of the CAR and AR scattering matrix elements. In the 
tunnel regime, the most probable process is normal re- 
flection. The EC scattering matrix element can be in- 
creased at the cost of the NR scattering matrix clement 
which has no signature in the conductance. But here, the 
AR contribution is as high as 3.3e 2 //i, while the maximal 
possible value in a four channel system is 4e 2 /h. 

If higher voltages are applied (see figure fT5b). each 
mountain range splits up into two less high mountain 
ranges. In the tunnel regime, the AR contribution disap- 
pears at higher voltage. Electrons and holes have differ- 
ent wavevectors and therefore different resonance lengths 
and if the resonances are narrow, electrons and holes are 
not in resonance at the same time. Here with the higher 
interface transparencies, the resonances are wider and if 
the electron is in resonance the hole part is still large 
enough to form a peak and vice versa. At higher applied 
voltage, the mountain ranges are surmounted by CAR 
peaks. The resonance of the EC elements of the scatter- 
ing matrix occur at different values of L r and CAR and 
EC are not in direct competition. The AR contribution 
(figure [T5] d) features at higher voltage two depressions, 
where there was only one at lower voltage. In one case 
AR is decreased in favor of CAR and in the second case, 
on the diagonals where L r = Li mod ir/kp, it is de- 
creased, as in the low energy case, in favor of EC. 




Figure 15: Differential conductance for symmetrically applied 
bias V = V a = Vb for intermediate interface transparency 
T = Ti nn = Ti ns = T rsn = T rnn = 0.5 and R = a) as 
a function of L r and Li at low voltage \e\V = 1 • 10~ 4 A, 

b) cut through a) along one of the mountain ranges (at 
Li — 3 183.77t/A;f) at higher resolution, c) as a function of L r 
and Li at higher voltage \e\V = 6.25 • 10 _2 A, d) cut through 

c) at Li — 3183.97r/fcF at higher resolution. 




Figure 16: Differential nonlocal conductance V a = for 
intermediate interface transparency T = 2] nn = T na = 
Trsn = Tmn = 0.5 and R = £, a) as a function of L r and 
Li at low voltage \e\Vb = 1 ■ 10~ 4 A, b) cut through a) at 
Li = 3 183.77r/fcF at higher resolution, c) as a function of L r 
and Li at higher voltage \e\Vb = 6.25 ■ 10 _2 A, d) cut through 
c) at Li — 3183.97r/fcF at higher resolution. 



2. Nonlocal conductance for Va = 0, Vb > 

Figure [TBI shows the nonlocal conductance. The peak 
positions at both studied voltages have the same symme- 
try as for the nonlocal conductance in the tunnel regime. 
But the peaks and depressions are less sharp and arc 
higher. In the tunnel regime, the EC and the CAR peak 
being superimposed at low voltage were similar in magni- 
tude. On figure [Tob however, the EC amplitude is much 
larger than the CAR amplitude. As in the symmetrically 
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Figure 17: Differential current cross-correlations for symmet- 
rically applied bias V = V a = VJ, for intermediate interface 
transparency T = Ti nn = T; ns = T rsn = T rnn = 0.5 and 

R = £, 

a) as a function of L r and Li at low voltage |e| V = 1-10 4 A, 

b) cut through a) along one of the mountain ranges (at 
Li = 3 183.77r/fc_F) at higher resolution, 

c) as a function of L r and Li at higher voltage 
\e\V = 6.25 • HT 2 A, 

d) cut through c) at Li = 3 183.57r/fcF at higher resolution, 
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Figure 18: Differential current cross-correlations for V a = 
for intermediate interface transparency T = Ti nn = T; ns = 
Trsn = Tmn = 0.5 and R = £, 

a) as a function of L r and Li at low voltage [e | Ve, = 1 ■ 10 _4 A, 

b) cut through a) at Li = 3 183.77t/A;f at higher resolution, 

c) cut through a) at L r = 3 183.77r/fcF at higher resolution, 
c) as a function of L r and Li at higher voltage | e[ Vj, = 6.25 • 
10 _2 A, 

dl) cut through c) at L; = 3 1 83.57r/fc_F at higher resolution, 
d2) cut through c) at L r = 3 183.77r/fcF at higher resolution. 



biased case, there is no CAR peak, but the forming peak 
features a crater. This figure illustrates that the dip in 
the CAR contribution occurs exactly at the peak position 
in the EC contribution, indicating that the competition 
between CAR and EC at the level of the scattering ma- 
trix is at the origin of the crater in the CAR peak. 



3. Differential current cross-correlations for a 
symmetrically applied bias 

Figure [17] shows the current cross-correlations for the 
symmetrical bias case. In the tunnel regime the cur- 
rent cross-correlations were dominated by CAR-NR and 
EC-NR as those processes depend only quadratically on 
the interface transparency, while MIXED has a cubic and 
CAR-AR, EC-AR even a bi-quadratic dependence on the 
interface transparency. For the current cross-correlations 
at intermediate interface transparency, this restriction 
does no longer exist. The current cross-correlations in fig- 
ure 117b feature trenches in L r and in Li direction. Their 
spatial extension lets us conclude that they are connected 
to the AR amplitude which is independent of L r respec- 
tively Li. From figure ITTb. we see that the trenches are 
due to MIXED. Figure [TTb shows also that the EC-AR 
component leads to sharp peaks at the positions where 
L r = Li mod ir/hp. These peaks are too narrow to be 
resolved in the surface plot [TTb ,. Here we have again an 
example of positive current cross-correlations which are 
not due to CAR. The CAR-NR component which was 
dominant in the tunnel regime is almost zero here, which 
is consistent with the small amplitude of the CAR pro- 
cess observed in the differential conductance. 

At higher energies (figures fTTH). the CAR-NR con- 
tribution gains importance as the scattering matrix el- 
ements for CAR and EC are no longer in resonance at 
the same position. In contrast to the tunnel regime, and 
like in the case of lower energies, the MIXED processes 
play an important role. 



4- Differential current cross-correlations for Va = 0, 
V b > 

The current cross-correlations in the nonlocal conduc- 
tance bias configuration shown in figure [18] feature at low 
energy, like the ones in the symmetrically biased case, 
trenches where the AR parts of the scattering matrix 
are in resonance. The MIXED component of the current 
cross-correlations plays an important role and leads to- 
gether with the EC-NR component to mainly negative 
current cross-correlations. At the points where L r = Li 
mod it/Icf and the EC elements of the scattering matrix 
are in resonance, the EC-AR component leads to positive 
peaks. At higher applied voltage, there is again the sepa- 
ration of the Li values where CAR or EC is in resonance, 
and MIXED becomes less important. 



14 



5. Discussion 

The above results show that it is not necessary to be in 
the tunneling limit in order to separate the different com- 
ponents of the conductance. Against the intuitive trend, 
the distance in the superconductor should not be too 
short. Owing to multiple scattering at the different in- 
terfaces, this can be understood as an optimum in terms 
of the CAR transmission, with respect to the normal bar- 
rier transmissions. In this case, the width and the height 
of the resonances in the conductance are wider than in 
the tunnel regime, which could even facilitate their obser- 
vation. As far as the current cross-correlations are con- 
cerned, they are richer but also more difficult to interpret 
in the regime of intermediate interface transparencies be- 
cause more components play a role. 



VI. COULD THESE EFFECTS BE OBSERVED 
EXPERIMENTALLY? 

What are the possibilities and obstacles for an exper- 
imental observation of the phenomena discussed in this 
chapter? The Fabry-Pcrot effect needs ballistic transport 
in a system with only a few (one) channels. The authors 
of the pioneering paper Ref. [30 observed Fabry- Perot in- 
terferences in single- walled nanotubes. Instead of varying 
the length of the nanotube, they change the Fermi level 
position in the nanotube by tuning the voltage applied 
to a gate under the nanotube. As the resonance posi- 
tion depends on the product of wavenumber and length, 
changing the Fermi wavelength has the same effect as 
changing the distance between barriers. This strategy 
is more realistic than to vary the distances Li and L r . 
However, if the electrons in the nanotube are sensitive to 
the electrostatic potential of a gate, they cannot be per- 
fectly screened. One would have to find a regime where 
the Coulomb interaction is small enough to be neglected 
in comparison to the, at least in the tunnel regime, also 
very small coupling energy between the nanotube and 
the reservoir, but high enough that the Fermi wave vec- 
tor can be influenced by a gate. 

One could imagine a device consisting of two normal 
conducting and a superconducting electrode connected 
by one long carbon nanotube with a symmetry similar to 
devices with nanowires used in recent experiment o 21 ' 3 5 ~ — 
where the normal and superconducting electrodes are 
connected over quantum dots. But while there, the sys- 
tem was operated in the Coulomb blockade regime, we 
need the other extreme of very small Coulomb energy. 
On both sides of the superconductor, there has to be a 
gate that can be separately tuned so that several peri- 
ods of Fabry-Perot interferences can be swept in order 
to obtain plots similar to the one presented in the last 
section. 

The next question concerns the interfaces. In Ref. l30l . 
where Fabry-Perot interferences are observed, the nan- 
otube is connected with near-perfect ohmic contacts to 



the electrodes. This implies a regime of high interface 
transparencies (T ss 1), not a tunnel regime. If the inter- 
face transparency is high, the requirement for the charg- 
ing energy to be much smaller than the coupling energy 
is easier to meet, but in the extreme case of T « 1 the 
contrast is too small. In the example of Refl30l the min- 
imal measured conductance is dl/dV ~ 3.le 2 /h and the 
maximal measured one dl/dV « 3.3e 2 //i. The last sec- 
tion showed that it is not necessary to be in the tun- 
nel regime to obtain a separation of the different con- 
tributions to the conductance. This can give hope that 
an intermediate value of interface transparency can be 
found. The next problem will then be to achieve this 
interface transparency experimentally. In addition, for 
sharp Fabry-Perot interferences, the length of the nan- 
otube has to be well defined. If one uses only one long 
nanotube the superconducting region is created by the 
proximity effect in the nanotube. So one cannot really 
speak of a sharp interface. 

In conclusion, it will be challenging to construct a de- 
vice where the separation of the different components of 
the conductance and of the current cross-correlations by 
interference effects in absence of electron-electron inter- 
action are observable. It could nevertheless be possible 
to reach an experimental regime where Fabry-Perot in- 
terference effects are in competition with Coulomb block- 
ade effects. The present analysis of a simplified model 
describing an extreme case could help to explain such 
experiments in an intermediate regime and it is a good 
basis for a more complete theory. 



VII. CONCLUSION 



We have shown that positive current cross-correlations 
in the tunnel regime cannot be enhanced with additional 
barriers by a process similar to reflcctionless tunneling, if 
an average (here, over the intcrbarrier lengths) has to be 
performed, mimicking disorder landscapes which are un- 
corrected on the two sides of the set-up. In analogy, one 
cannot therefore expect that the use of diffusive normal 
metals will facilitate the experimental observation of pos- 
itive current cross-correlations. This negative conclusion 
is important for experiments. 

On the contrary, the insertion of additional tunnel bar- 
riers into ballistic NSN-devices leads to Fabry-Perot like 
oscillations of the differential conductance and of the cur- 
rent cross-correlations. The resonances for different pro- 
cesses are well separated at finite energy, because of the 
different energy dependence of electron and hole wave 
vectors and CAR and EC processes can be efficiently fil- 
tered for intermediate transparencies and for a central su- 
perconductor width of the order of the coherence length. 
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Appendix A: Details on the scattering approach 

The elements s*^ of the scattering matrix are cal- 
culated within the BTK approach^: Two-component 
wave functions, where the upper component describes 
electrons and the lower components holes, are matched at 
the interfaces and the coefficients of the resulting system 
of equations give the elements of the scattering matrix. 
In the normal conductors the the wavefunctions are plane 
waves with wavevectors close to the Fermi wavevector 
Kk F = y/2mfi. The wavevector for electrons reads Hq + = 
V2m^/fi + E, the one for holes hq + = s/2rn^/Jl — E. In 
the superconductor, the wavefunction has to obey the 



ducting gap is supposed to be a positive constant in- 
side the superconductor and zero outside of it. This is 
achieved by modifying the amplitudes of the wavefunc- 
tion in the superconductor with the coherence factors 
ue and ve which read for energies smaller than the gap 
(E < A): 



WA 2 - E 2 
E ' 



VE 



7T 



WA 2 - E 2 
E 

(Al) 



1 - 



and by using for quasi-particlcs proportional to 



the wavevector hk + = \J1m\/ [i + iy/ A 2 — E 2 and for 



quasi-particles proportional to 



ue 



the wavevector 



Hk~ = \[ T Zm\/ fi — W A 2 — E 2 . 

For example, the wavefunctions for an electron incoming 



Bogoliubov-De Gennes equation, where the supercon- from electrode N a take the form 
I 



1 e lq x + sf a e~ lq x 
ci e lq+x +c 2 e 



c 3 e lq x + c 4 e - iq ~ (x - Ll) 



(?) (*aa e«-* + e- iq ~*) 

tq+{X ~ Ll) ) + (!) ( 
(S) ( C5 e%k+[X ~ Ll) + c e e- ik+ ^~ L '- R ^ 

(^) ( Cy e -«-(*-L, +Cg e *fc-(x-L I -fl)^ ^ 

( C9 e ** + (*-L>-m +Cl0 e -il + {x- Ll -R-L r )\ 

(c u e icr{x - Ll - R) +ci 2 e -^(^-L,-R-L 

) 



iq~(x~L,-R) + <A e e -iq-(x-L t -R-L r y 



(A2) 

(A3) 
(A4) 
(A5) 
(A6) 



in the sections N a , N/, S, N r , and N& respectively [see 
Fig. Q] and give access to the scattering matrix elements 

s a a j s a,a: s b'a an d s b'a- The remaining elements of the 
scattering matrix can be obtained from the other possible 
scattering processes i.e. a hole incoming from electrode 
N a , an electron/hole incoming from electrode Nf,. 

The interfaces are modeled by <5-potentials 
V{x) = Zhv F S(x), where the BTK parameter Z 
is connected to the interface transparency T by 
T = (1 + Z 2 )" 1 . The elements of the scattering matrix 
can be determined and the constants c, eliminated using 



the continuity of the wavefunctions at the interfaces 
[^JVa(O) = V'JVi(O) etc.] and the boundary condition for 
the derivatives [V'jv^O) — fp' Na {0) = Zhv F ^pa(0) etc.] at 
every interface. 

In simple cases, i. e. for only two or three sections 
and in the limit of zero energy, the system of equations 
giving the scattering matrix elements can be solved an- 
alytically (seei£), but in the present case of five section 
the expressions become so unhandy that we decided to 
solve the the equations numerically 
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Appendix B: Components of Current Cross 
Correlations 



Components of the current cross-correlations: 

o,v a ,v b) = 2 -^f dE ( 



S ab< T 
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Scattering matrix 
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Green's function 
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CAR-NR 


CAR 


EC-AR 


AR-AR 


MIXED1, MIXED 2 


PRIME 


EC-NR 


EC 


CAR-AR 


AR-AR 


MIXED3, MIXED4 


MIXED 



Table I: Correspondences between the categories in the lan- 
guage of Green's functions from 1 ^ and in the language of scat- 
tering matrix elements. 



Appendix C: Relations between the noise 
classification in the BTK and in the Green's 
functions approach 

The elements of the scattering-matrix are connected to 
the retarded Green's functions of the tight binding model 



studied h>!£ via 



8?f = My + 2^ x 7:\ i>:< (ci) 



where U is the transmission coefficient of the barrier i, 
pf the density of electron or hole states of electrode i 
and Gija/s the Green's function connecting the first site 
in the superconductor next to the electrode j to the first 
site in the superconductor next to the electrode i. 
Table|T]shows the correspondences between the categories 
in the language of Green's functions and in the language 
of scattering matrix elements. 
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